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ABSTRACT 


The rjyy model for the equatorial electrojet presented earlier is 
revised by including a meridional currant system, thus making the currant 
density divorgence-froe* Using a spherical harmonic expansion for the 
geomagnetic field a model electro jet is constructed for each of the longi- 
tude h of n*% /| 0 °j g(} 0 | igo^j and 280°R, It is shown that the en .t-ntial 
feature »#j Uus murid tonal current system is the existence o f two current 
lonjm, **,» on each side of the dip equator and with an cquatorward currant 
in tin* lower K region, an upward current over the dip equator, and a re- 
turn current flowing in the upper E region. The intensities of t m 
electroji't and the associated meridional current are significantly 
depend nt on longitude, being strongest over Peru and weakest over India, 
The new model predicts a reversed (westward) current in a thin, narrow 
layer located at several degrees off the dip equator, which is due to a 
rever»al of the vertical polarization field. The meridional current 
loop creates a toroidal magnetic field of an intensity comparable with 
that of the main electro jot. 



INTRODUCTION 


« 


The equatorial electro jet has been a subject of extensive study both 
observational and theoretical* Direct observation* *f its magnetic field 
by rocket-borm* magnetometers have been made by Singer et al * (1951) and 
Cali ill Cl 959) in the Pacific, by Maynard et al. (1965), Maynard and 
Cahill (1965), and Saatry (1968) over Thumba, India, and by Davis et al* 
ft967) and Maynard (1967) off the coast of Peru* Numerous papers have 
been written on ground observations of the electro jet field j and a summery 
with extensive references has been given by Onwume chilli (1967), 

Sugiura and Cain (1966) presented a model lor the equatorial electro** 
jet which was based on the assumption that the vertical Hall current is 
completely inhibited everywhere by the polarisation electric field. With 


this assumption the electrojot current density is proportional to the 


effective conductivity ftyy, 


in the east-west direction, if the driving 


electrostatic field is uniform and eastward. This model may be regarded 


as the xero-th order approximation, since the current density is simply 
calculated using a set of given variables without solving the differential 
equation that a self-consistent model must satisfy, A merit of the approx- 


imation is that the realistic magnetic field configuration can readily be 
incorporated. Indeed, the main purpose of their paper, which will be 
referred to below as paper 1, was to investigate the longitudinal variation 
of the equatorial olectrojeu by using the magnetic field calculated from 
the spherical harmonic expansion of Jensen and Cain (1962) j the height 
distributions of the electron density and other relevant atmospheric para- 
meters were assumed to be independent of longitude. 
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However | being a ziro-th order approximation the model doe* not 
exactly satisfy the differential equation to be discussed below 5 stated 
more explicitly, the current density, j, is not divergence- free, i,e M 
^*j 5^0. Untiedt (1967) has obtained a solution that satisfies the con- 
dition 9 «J® 0 j his solution includes vertical currants and consequently a 
meridional current system. The model, however, assumes a dipole magnetic 
field. 

The present paper gives an Improved model which satisfies v» j*0 and 

it 

which at the same time usee, though in a regional sense, the magnetic 
field configuration represented by a spherical harmonic expansion. As a 
new feature the model indicates that there can be a reversed, westward 
current in a thin layer at several degrees off, and on each side of, the 
dip equator. The reversed currents are essentially due to a reversal of 
the vortical polarization field rather than being a ’return current 1 from 
the main eastward electro jet. The magnetic field measurement made by 
Maynard (1967) with a magnetometer flown on a rocket off the coast of Peru 
has shown the existence of a reversed current similar to that predicted 
in the present model, 

THE COORDINATE SYSTEM AND THE ASSUMPTIONS 

Geocentric spherical coordinates (r, 0, d) are used in this paper 
for the convenience of calculating the magnetic field from the spherical 
harmonic expansion given by Cain et al, (1967)? the field is calculated 
for epoch 1965.0. The spherical coordinate system is defined in the fol- 
lowing manner. Let the geographic longitude at which we wish to obtain 
a model for the equatorial electrojet be First, the geographic 
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colftFUude. 1 i, of tin? dip equator at this longitude and the direction of 
the tangent to the dip equator at the position Aj), both on tho 
earth *n mirfato, art* determined from th§ apharietl harmonic expansion of 
the magnetic field, Then a (geocentric) spherical coordinate system in 
defined an that the tangent to the dip equator lies in it a equatorial 
plane (Q*¥t/2)s spherical coordinate® (r, 0, 4) used below refer to thin 
ays tom* Let the longitude of the point (%, A|) be 4\ in the new spherical 
coordinate system, The crocs section of the electrojet and the meridional 
current system are determined in the meridian plane of 4*4 The magnitude 
and the direction of the magnetic field are calculated from the spherical 
harmonic expansion in this meridian plane ai function* of r and 0, Once 
the magnetic field is calculated in this meridian plane we assume that 
the magnetic field in the model calculation in Independent of 4* The same 
assumption is made with ail other quantities chat determine the conductivity. 
It i& assumed that the magnetic field vector lies in the meridian plane 
Ml this is a justifiable assumption near the dip equator. Since the 
model is Intended to represent the equatorial electrojet regionally, a 
two-dimensional model is considered to be satisfactory. The model elec- 
trojet is assumed to be in a steady state, 

TUB FORMULATION 

The basic equations that the electric field intensity, E, and the 
current density, j, satisfy are 

y x E ■ o 
7 • j “ 0 


( 1 ) 

( 2 ) 


Equation Cl) follows from the trendy state condition* and (2) expresses 
Che divergence*' free current, With the condition Cl) gives ¥,$ 

mil m equation lor E r and Egi 


and 


Tc (rE 0> 



(&) 


where E^ i§ a constant, and a in the radius ol the earth. Within the 
boundaries considered hare is nearly constant. In Untledt f » recti- 
linear (plane earth) model this correspond# to the constant Ey , 

The components of j can be written as 

Jr ■ + «,co» n I)E E + (c 0 - *5,) Ee*inIco*I - <?*8^eo»I <5) 

jp - (n 0 - u, ) E r ,lnIco»X + (sr 0 co» M I + ■7j«ln : 'I) E ft + sr t ,E,J»inI (6) 

Jd " 'Tsfirco*! " «f->Efl*inI + a, E,j (7) 

where a 0 » , and <jr s are the direct, Pedersen, and Hall conductivities, 

respectively, and 1 is inclination of the magnetic field. Equation (7) 
is to be used to calculate when and E@ are determined from (2), (3), 
(4), (5), and (6). 

Since 0, it follows from (2) that in the r-0 plane and Jq 

can be derived from a current function f such that 



fttm tSi and th) with tl) t t f rM ft* nin b** aspr*****! in t*»r*«* 

oi ! r anti j 4 % a i| hence a* I tint Hons of derivative* of f* &«h*t Util in? 
tlKMM* in th% m obtain i different tal eqimtitift for f a* follows* 


r* £-X * 24 ;££_ + *, zjl + l 

«r »r in #>n* % 


? r . K * f* * O HO) 


whero f| » if t *,»*, fn ar > function* of r and © and are given by 


i i * Po cor 1 ! * Oj a in' I 


fg » r wi fM| * p 0 ) *ln i oo* I 
f . * r** f> 0 sin* r * §i co* fe I) 

f 4 • r* 1 ^ % cor I *1* ^airft) * y tCbj - bo) sinl eosll 

• Cm*- p t> ) ainl coal cot©] 


k * r"*[r ^ { (bj -po) *inl coil] + (v 0 a in' I *l* pj »;or 1) 

- Cbj-Mo) »fnl coal * <p 0 sin ; 'X A* cor*I) cot 9] 

U - - r* 1 lr ^ <2t iln i) + ^ (S, cod) - £* eo*l eotfl) Rfo 
where * U®o and p t «* l/*?| * 


THE BOUNDARY CONDITIONS 

Wo impose the following boundary condition* t (1) f»0 at the lower 
(Wi> and the upper (r»*r 8 ) boundaries, (ii) At the northern (@*0^) and 
the southern (9-n**9 1 ) boundaries, which are far enough from the dip e qua to 
j r “Q and the ©-component of the height-integrated current is equal to the 


current that would be expected if the conducting layer were infinitely 


thin and having a conductivity equal to the height-integrated conductivity 



of tin? lay* r, Wwe tin* drivinn electric f l« IU Um o «»ily *r. 

in the cast* treated here, kUo latter condition reduce** to shut the &*tm* 
pommt of tin hoighc*into£r«6ed current vanishes, It in noted that Eli# 
above boundary conditions hero** the sssmi m those used by IMiedt CiW) 
in th# plane-earth limit. 

The condition (if) enables m to express Eg in forms of K$ 0 , and 
f at d*0| and n * © 4 becomes 





dr 


%iin I 4* tyvoaPi 

fo gp^iinX 
4| i * Ojcpi^ 



rg 00^1 
^ hoSi^ f * O t COS* I 


X 


V>1 


Jr t e 0 *in T I + n t eor'X 


rdr 


€12 


Equation (12) shows that f vanishes at r»r x , but that % Is not exactly 
sera at r*r^ . However, f is nearly xaro at r*r PJ> because the two terns 
in (12) approximately cancel each other, and, in addition, because ty% md 
ff % are both nearly sero at r“r 4 and r s , We could make f exactly aero at 
r*rg by making o x and % vanish at this height, T.t practice, however, 
there is little difference between the two treatments. 

In the actual calculation r x and r§ are chosen to represent the 
heights of 70 km and 380 km, respectively, and 9 % to be 80°, or 10° off 
the dip equator. 
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t m $mntm 

Fhv di»t-r tfetnent l|*f| f* tn CIO) l* found to bv negat ive, and ihv 
*Uf r»r***ul«i ©quel i*m i* a? elliptic type, Equation (10) ha* been solved 
nxmvti* 4t iy toy »*i«K thr method of successive ovsr-relaxat ion* tn r*«* 
film In* Mu) hy a ** fc ot different* equations, grid point* were spaced ut 
jot* rv4i* »»? % km t«* r and 0.5° in 0* The partial derivatives of given 

r * in iht •urn icintt* f 4 , fa, f i were calculated nua&r leal iy using 
4 f iiirr «< *h« 

Hit* value iitr T.ja was taken to be 2,4 x 10*^ volts/®, which ut the 
»*»© 4* in paper ! . the solution i* propcn tional to E^ot hence it 
t en tot* rend! ly adjusted to different values of £^ 0 . The electron density 
profile end ilw height distribution* of o Q , o % , end ere given in the 
Ap p^otUx * Solutions were obtained for longitude* 0°, 40°, 80^, 180^, and 
280°p , roughly representing Eaet and West Africa, India, tne South Pacific, 
and Smith America, respectively * Figure 1 show* contour* of equal t for 
these longitudes* Contour* are drawn so that the current flowing In a 

(meridional) slab of thickness 1 km and between successive contour* 1* 

* 

}Q amperes, Th© direction of the current I* upward between the two foci* 

F l ipiro 1 indicates a strong longitude dependence of the intensity of 
the mertuioniti current, To make a quantitative comparison the total current 
(per mei^r- thickness) flowing between the two foci is given in Table 1 for 
live different longitudes; in preparing Table 1 the driving eastward electric 
field, K^ (> , i* normalised to 1.0 x 10“ 3 volts/ra, 1 * 0 *, 1/2,4 time* the 
value used for Figure 1, It can be teen fro® Table 1 that the meridional 
current i* strongest over South America and weakest over India; this re* 

1 at Ion is opposite to that of the strength of the earth 1 * magnetic field, namely, 
the latter is weakest over South , ‘merica and strongest over India, 



Tabi* I. Th# total current (in amparaa par mts r» 
ehickno*«) flowing batwaan tha two foci and tha rati© 
to tlia total currant for 280®&t » 1.0 x 10”^ volfca/m. 


Geographic jpngUuda 

Total currant 

Hatio 

0° K 

0.224 

0.76 

*0° 

0.191 

0.65 

80° 

0.166 

0.56 

180° 

0.228 

0.77 

280° 

0.296 

X.00 


m 
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THE BLECTROJET 

By numerically differentiating ¥, \ v and Ja can h® calculated from 
(8) and (9); and the meridional component* of electric fit* Id tan be 
obtained by solving (5) and (8) for E r and Bn* Then the eastward current 
Jd can be determined from (7), Figure 2 shows contour* of equal current 
intensity (in units of 10“^ amps/m^) for longitude 280°E ; the driving 
electric field E^ 0 in taken to be 2,4 x 10* * volts/m, The peak current 
density is 3*4 x 10*^ amps/m^, and In located at about 105 km altitude 
and approximately 0 * 1° to the south of th<s surface dip equator* The 
center of the electro jet approximately coincides with the dip equator as 
defined at tl\e height of the center; at the height of 110 km the dip 
equator is approximately 0*3° south of that at the ground for this longi- 
tude. 

* 

An unexpected feature in the cross section is the existence of 
two thin, narrow regions at several degrees off, and on both sides of, 
the dip equator, in which the direction of the current is reversed, To 
ascertain that these are not from effects from the boundary conditions 
at Latitudes ± 10°, the whole calculation Was repeated with a new boundary 
condition that ¥*0 at these boundaries. Not only did we obtain a meridional 
current system that is nearly identical except near the boundaries where 
the current is weak, we confirmed that the essential features of the j^ 
profile is the same as in Figure 2, The region of the westward current 
to the south of the electro jet is larger and its intensity is stronger 
than the northern counterpart, This must be a result of the asymmetry 
of the magnetic field at this longitude, The maximum westward current 
intensity, located at 109 km altitude and 7° south of the dip squalor, 


/ 


ie 0*13 x 10"=* ampa/m^, and 1 b very much weaker than the maximum electro- 
jet current density. Nevertheless the existence of the reversed currents 
is of interest, because Maynard (1967) has shown a rocket measurement of 
the magnetic field that suggests the presence of a reversed current at a 
location several degrees south of the dip equator off the coast of Peru. 

The cause for the reversal in the currant lies in the reversal of the 
vertical polarisation field E r , as is expected from (7)» Figure 3 shows 
the height distributions of E r for longitude 280°E . At 5° north and south 
of the dip equator, E r has a very steep gradient at the altitudes of about 
104 km and 114 km resulting in a hole, and a reversal, in the polarization 
field. The depth of the hole Is found to depend on longitude} for the 
longitude of 80°B, the polarization fiald decreases roughly in the same 
region, but the decrease is not large enough to invert the field except 
in a very small area} and consequently is likewise not reversed, only 
having a steep depression, except in the small area of the inverted electric 
field. Figures 4 and 5 show the cross section of and the vertical polar- 
ization field, E r , for longitude 80°E, Comparing Figure 4 with Figure 2 
it is seen that the electro jet over India is expected to be weaker and 
its width tc be narrower than over South America, 

TOROIDAL MAGNETIC FIELDS 

In a steady state, B is related to J by the equation 

V x B * 4rrj 

From this relation together with h/bd * 0 we obtain 


(r sin 0/4TTa) B^ + “ constant 










'Hie value of f i* zero afe r*t % and and ii nearly zero afc the boundaries 
0»% and 11 *% . Hence m can determine approximate intensity of the toroidal 
magnetic field relative to B^ at the boundary. Since r/a « l and iin Q n l, 
approximately 

n w An»|f *1 B^o (16) 

where B^ 0 is Bg{ at the boundary. 

<*1 

If the ^ value* indicated In Figure 1 are multiplied by 4t* x 10 ,> 1.25* 

in y is obtained, The field ii westward to the north of the dip equator 
and is eastward to Hie south. The maximum value* of relative to the 
boundary value, for the longitude* of 0°, 40°, 80°, 180°, and 280°E are 
given in Table 2, where the Intensities are normalized to a driving 
electric field of 1,0 x IQ"** volta/m as in Table 1. The magnitude of the 
toroidal magnetic field at the center of each current loop is comparable 
with that of the electro jet field as observed on the ground. It will be 
of interest to directly observe the toroidal fields by a vector field 
measurement with a rocket-borne instrument. Rocket observations so far 
made have all been scalar measurements, 

DISCUSSIONS 

(1) The meridional current system presented in this paper differs 
front that given by Untiedt (1967) in that in our system a current loop is 
formed on each side of the dip equator, while in Untiedt’ s model the 
meridional current flows into the equators! region from higher latitudes 
at the lower dynamo altitudes and flows out again at higher altitudes: 
see his Figure 2, 

The difference between Untledt’s and our results is thought to stem 
mainly from the following circumstance. Assuming, for simplicity, a dipole 
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field configuration, the tine of magnetic Force that i *m lit* * the «|»iht 
boundary in our model (i*e,, 380 km height) over tin *iip equator hwinivj 
the boundary nt about 12° latitude, which in very etoae to tin* nor t hern 
(or southern) boundary in our calculation* While* in Ontiedt: *n model the 
line of force touching hit upper boundary (300 km height) over the equator 
is cut off by the vertical boundary at about half way to its intersection 
with the lower boundary* Since the conductivity a 0 along lines of force 
is very large, we believe that the width of 500 km in Untiedt*# model is 
not quite adequate to obtain a reliable meridional current ayatem. how- 
ever, hl« model doss agree with ours in gross aspect# in the region with- 
in 2° or 3° of the dip equator* 

(2) To compare the distribution in the new model with that in 

the Oyy model the cross section of Oyy ia plotted in Figure 6 for the 

longitude of 280°E * The half-width, as defined by the distance from the 

peak to one-half the peak, is about 0.7° in the o V y model, while the 

corresponding half-width in the new model is approximately 2,5°. The 

latter value (** 250 km) is in rough agreement with an estimate from Figure 

5 of Untiedt’s paper. Over India the half-width la estimated to be about 

**0 

2°. In Xigure 6 the peak value of cjyy is 1.97 x 10 mho/m. With the 
driving electric field of 2.4 x 10" 3 volts/m the maximum current density 
in the Oyy model is 4.7 x 10“ 5 amps/m 2 , which is slightly less than the 
peak current density of 5.4 x 10"^ amps/m 2 in the new model (Figure 2). 



ia 


(i> Referring to Figure $ eh* maximum vertical polarisation Halt! 
over the dip equator i* at £ht height of 100 km and it* value ia approxi- 
mately 25 time* the driving eastward field* This large ratio of the 
vertical field to the eastward field ia the ««aene© of the existence of 
the equatorial electro jet * In the Oyy model this ratio ia n%ln % over 
the dip equator | and the ratio is o maximum at 100 km height, its value 
being 28, Hence the meridional current somewhat reduces the vertical 
♦•It i trie field hut only by it small amount, 

f4) In the model presented in this paper the driving electric field 
is assumed to be in the east -wait direction. Effects of an addition of 


a north-south component of the electric field were investigated by including 
this component in the formulation, hut the results indicated that this 
doas not alter the essential features of the meridional current system, 

A similar conclusion has been reached by Untiedt (1967) in his model 
calculation. 


(5) The maximum electrojet current density given here is about 
5 x 1CT 5 ampa/m 2 for the longitude of 280°E* The corresponding quantity 
in Untiedt f s model is about 0*9 x 10**5 amps/m 2 * The driving electric field 


used in the present model is 2,4 times the value assumed by Untiedt , The 
difference in the ratio a%h x in the two models may account for the re- 


maining factor of about 2, though this point must be examined further in 
the future* When the rocket measurements were made off the coast of 
Peru by Maynard (1967) and Davis et al* (1967), the maximum current density 
over the dip equator was approximately 0.7 x IQ" 5 amps/m 2 * Sastry (1968) 
measured a maximum current density of approximately 4 x 10**^ amps/m^ over 
Thumbs, India, Thus in both cases the observed values are smaller than 


the theoretical value* given here for the rcapectivc longifcydt** hy the 
§&m factor of about ?* Thii difference ii thought to be dm* mainly to 
weaker electric field* prevailing at the time* of the rocket measurement* 


than aaiumed in the model, but difference* in the electron d*n*Hy prof i l 
and other ionospheric parameter* may contribute to a certain extent * 
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. the How pattern of the meridional current anauciated with 

the equatorial electro jet for different longitudes, the 4 urnnt 
flowing between successive contours is 50 amperes per krr hick- 
mess. The driving electric field, E<| 0 »2,4 x 10“ ^ volts/m. The 
current flows upward over the dip equator. 

The eastward current density, for longitude 280°E ; In units 
of 10"* amp/m 2 , The driving electric field, E^ 0 * 2,4 x 10*3 

: 

vtilts/ffl. The shaded areas represent reversed, westward currents, 

, The vertical polarisation field, Ej., directly over the dip 
equator and at 5° north and south of the dip equator; for 2800 f 
longitude, flu* driving electric field, K^y* 2,4 x 10“^ vol ts/m* 

, i in 1 . astwurd current dunaity, for longitude 8 Q°Et in units 
o f 10“ * amp/m 2 . The driving electric field, E^ y *■ 2*4 x 10"^ 
volta/m. The shaded areas represent reversed, westward currents, 

. The vertical polarisation field, directly over the dip 

1 

equator and at 5° north and south of the dip equator t for 80°E j 

longitude; E ^ 0 * 2,4 x 10“ 2 volts /m, J 

J 

, Tin* ^yy profile; in units of 10“^ mho/m. The peak value is 1 

1,97 x 10“ 2 mho/m, j 

, The electron density profile used in the model, I 

. Conductivities o 0 , and ^2 aa functions of height, over the j 

dip equator and for the longitude of 280°E, j 
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